Non-topographic description of inherent structure dynamics in glass formers 
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We show that the dynamics between inherent structures in glass forming systems can be under- 
stood in purely dynamical terms, without any reference to "topographic" features of the potential 
energy landscape. This "non-topographic" interpretation is based instead on the existence of dy- 
namical heterogeneities and on their statistical properties. Our view is supported by the study of 
simple dynamically facilitated models of glass formers. These models also allow for the formulation 
of quantitative theoretical predictions which are successfully compared to published data obtained 
in numerical and experimental studies of local dynamics of supercooled liquids. 
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It's been like a kind of capsule, a bubble in time and space... 

D. Lodge, Paradise News 



I. INTRODUCTION 

In 1969, Goldstein suggested that the dynamical slow- 
down of a liquid approaching its glass transition could 
be understood in terms of its potential energy land- 
scape pj . The practical implementation of this descrip- 
tion was later introduced by Stillinger and Weber via 
the inherent structure (IS) formalism ,2|. The landscape 
or "topographic" view has become one of the paradigms 
in the field 0] . The purpose of this paper is to pur- 
sue, however, the rather orthogonal idea that the land- 
scape description is essentially incomplete, in the sense 
that from the statistical properties of IS, or similar static 
characteristics, it is not possible to capture the central 
spatial and dynamical aspects of the physics of glass for- 
mers. Instead, we aim here at a "non-topographic" rein- 
terpretation of the dynamics between IS, showing that 
the central role is played by dynamical heterogeneities, 
without resorting to any special features of the potential 
energy landscape. 

The connection between the existence of a dynamical 
correlation length and the landscape description is in fact 
implicitly present in the original papers 0,0, althou gh i t 
only clearly appears in more recent works [a, S H l&llOfl • 
The connection relies on two remarks, (i) According to 
Goldstein Q], the landscape is useful for "viscous" liq- 
uids, for which a separation of time scales between fast 
vibrations and slower structural relaxations exists. This 
gives rise to the notions of "basins" or "traps" which are 
then purely dynamical concepts, (ii) Escape from these 
traps was described to be "localized" in space pj, a fact 
which was later numerically confirmed [5|- The tran- 
sition rate between basins therefore becomes extensive 
"for large enough systems" p|. How large the system 
must be, what is the geometry, and what are the typ- 
ical length and time scales of these rearrangements are 
crucial questions left unanswered unless spatial and dy- 
namical aspects are included explicitly in the theoretical 



description. 

In what follows we show that natural answers can be 
obtained in the context of dynamical heterogeneities i 
thus leading to a consistent qualitative interpretation of 
IS dynamics. This we demonstrate using the simplest 
models which capture the phenomenon of dynamic het- 
erogeneity, and for which the landscape properties are 
completely irrelevant (sections II and III). We also show 
that a satisfactory quantitative agreement with numeri- 
cal and experimental studies of mesoscopic systems can 
be obtained (section IV), emphasizing the relevance of 
the non-topographic description of supercooled liquid dy- 
namics discussed in this paper. 



II. DYNAMICAL HETEROGENEITIES VERSUS 
"HIDDEN INHERENT STRUCTURES" 

In a cold dense liquid mobility is sparse. A microscopic 
immobile region is mostly surrounded by other immobile 
regions. It can become mobile only if it is next to a rare 
microscopic mobile region: its dynamics is facilitated by 
the presence of an unjammed mobile neighbour. This is 
the concept of dynamic facilitation, originally introduced 
in [T^ | , which leads to the phenomenon of dynamic het- 
erogeneity 0,0]. 

The elementary models which capture the idea of dy- 
namic facilitation are the Fredrickson- Andersen (FA) [lj 
and the East 0] spin facilitated models (see ^(| for a 
comprehensive review). They consist of a chain of two- 
state spins, rii = 0, 1 (i = 1 • • ■ N), with non-interacting 
Hamiltonian H — an d single spin flip dynamics 

subject to local kinetic constraints. In the FA model, a 
spin can flip if either of its nearest neighbours is in the 
up state, while in the East model a spin can flip only 
if its nearest neighbour to the right is up. The equilib- 
rium behaviour of both models is that of an ideal gas of 
binary excitations, with an equilibrium concentration of 
up spins given by c = 1/ (l + e 1 ^ T ^ at temperature T. At 
low temperatures the dynamics is glassy due to the com- 
petition between decreasing the energy and the need for 
facilitating spins. The FA model is Arrhenius, or strong, 
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FIG. 1: Top: an IS trajectory in the FA model at T = 0.6 
and system size L — 200 = 32£(T); vertical axis is space, 
horizontal one is time, up spins are black. Bottom: the cor- 
responding time series of the IS energy. 



in the jargon of supercooled liquids [Fjjj with relaxation 
time t rc i = e 3 / T . The East model is super- Arrhenius, or 
fragile, t Tei = e 1 '^ 2 ln2 l 

In order to study the dynamics of IS in the FA and East 
models we consider the single spin Monte Carlo dynam- 
ics in equilibrium at temperature T, using a combination 
of Metropolis and continuous time algorithms 0]. At 
each Monte Carlo step, we quench the system by run- 
ning a zero temperature Metropolis dynamics. We thus 
obtain the finite T and its corresponding IS trajectory, 
as initially suggested in the context of liquids. In Fig. 
(top), we show such an IS trajectory in the FA model at 
T = 0.6. The vertical axis corresponds to space and the 
horizontal one to time. Up (down) spins in the IS con- 
figuration are denoted black (white). The IS structure 
trajectory is very close to that of the corresponding fi- 
nite T one. This is due to the fact that, as a consequence 
of the kinetic constraints, an up spin can be relaxed un- 
der the quench only if a neighbour is in the up state, so 
that the quenching procedure does not break the con- 
tinuity of the up spin "world-lines". This IS trajectory 
then displays the structure of a dense mixture of "bub- 
bles" separated by excitation lines, as discussed in [l4| . 
In Fig. HI (bottom), we show the evolution of eis, the 
(extensive) energy of the IS, as a function of time t. The 
time series is featureless, corresponding to a sequence of 
many uncorrelated events. This is analogous to what is 
found in simulations of liquids for high temperatures or 




FIG. 2: Same as Fig. 1 but for T = 0.4 and L = 53 = 4l(T). 



large system sizes 0, El • 

The existence of a dynamical coherence length, £(T), 
is however clear from the IS trajectory in Fig.^ Individ- 
ual transitions in IS correspond to the birth or closure 
of the trajectory space bubbles, i.e., the branching or co- 
alescence of excitation lines. The dynamical correlation 
length £(T) is thus given by the typical size of the bub- 
bles, which is fixed by the equilibrium concentration of 
excitations, £(T) = c _1 (T). In Fig. ^ the system size 
L is much larger than £(T), L = 32£(T). In Fig. we 
show instead an IS trajectory and its corresponding e/s 
time series in a case where the system size is much closer 
to the coherence length, L = i£(T). The time series for 
eis now displays the coherent changes in IS energies also 
observed in numerical simulations of "sufficiently small" 
supercooled liquids la, Uj|- Obviously, "sufficiently 
small" only becomes a well-defined statement when £(T) 
is first defined. We emphasize also that these "hidden 
structures" follow here from the presence of dynamical 
heterogeneities, with no obvious relation to specific to- 
pographic features of the energy landscape, which in this 
case is trivial. 

All the features of the e/s time series in Fig. [21 (bottom) 
can be traced back to the slow bubbles, i.e., dynamic het- 
erogeneities, in the trajectory of Fig. [21 (top). The fast 
spikes in Fig. [21 (bottom) correspond to smaller bubbles 
that "wet" the larger ones, which in turn determine the 
more persistent plateaus in e/<j. In principle one could 
renormalize out the smaller bubbles and obtain a time 
series for "metabasins" in analogy with H, EJ • In our 
case this is unnecessar y s ince we know the space time 
distribution of bubbles [lj| and we can therefore account 
for the effect in the IS dynamics of the whole range of 
dynamical heterogeneities. This is especially important 
in the case of the East model, where the hierarchical na- 
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ture of the dynamics leads to a fractal wetting structure 
which in turn is responsible for the temperature depen- 
dent stretching of correlation functions characteristic of 
fragile systems 0, |2(J • 



III. LIFETIME OF HETEROGENEITIES 
VERSUS "HOPPING TIMES" 

Each event in the time series of e/s of Figs. and 
121 which increases/decreases the energy of the IS cor- 
responds to the birth/closure of a bubble in trajectory 
space. The distribution of time intervals between events 
is then given by the distribution of time extensions of 
bubbles, or lifetimes of dynamical heterogeneities, p(t). 
One has [3: 



P(*) 



p(£)p(t\£)d£ : 



(i) 



where p(£) — ce~ ci is the distribution of length scales of 
heterogeneities, p(t\£) the conditional distribution of cor- 
responding time scales, and J\f a normalization factor. 
The stretching exponent (3 is temperature independent, 
(3 = 1/2, for the FA model, and decreases with T for the 
East model, with (3 cx T at low temperatures. The corre- 
lation functions are determined through the persistence 
function P(t) = J t °°p{t')dt' ~ exp[-(i/iroi) /3 ]- These 
models are an explicit realization of stretched relaxation 
induced by the distribution of time and length scales of 
heterogeneous regions 

To numerically access the distribution QJ, we follow 
the empirical procedure of Refs. [|| Q . As for these sim- 
ulations of supercooled liquids, we strike a compromise 
between small enough system sizes, comparable to the 
dynamic correlation length in order to distinguish be- 
tween individual events, with large enough ones to avoid 
finite size effects, corresponding here to an upper cutoff 
of the integral JQ). Again, the procedure becomes well- 
defined when p(£) has first been defined, Eq. J[J. In prac- 
tical terms, we seek the smallest system size for which 
autocorrelation functions converge to those of the bulk. 
In the case of the FA and East models we have checked 
that this is obtained for systems of length L ~ 4£(T) 
and 8£(T), respectively. The difference stems from the 
hierarchical structure of time scales of the East model 
which enhances the relevance of larger length scales with 
respect to the FA case. 

Figs. and Q] show that Eq. indeed describes the 
distribution of times between events — or "hopping times" 
in the language of H, || — in the IS time series. Fig. |3| 
(top) presents for the FA model, at various low tem- 
peratures for system sizes L — A£(T), the IS energy 
correlation function, C(t) = [(eis(t)eis(O)) — (ej s )]/L. 
An almost temperature independent stretching, close to 
j3 = 1/2, is obtained, as expected at low temperatures, 
Eq. QJ. In the lower panel of Fig. [3] we show the dis- 
tribution of the logarithm of hopping times, 7r(lnt), for 





FIG. 3: Top: IS energy correlation function C(t) for the FA 
model at various T = 1.0, 0.6, 0.4, 0.3 and 0.2 (from left 
to right) and system sizes L = 4£(T), in a log-double log 
scale. The lines indicate stretching exponents /3 — 0.57 and 
(3 = 0.5 (from left to right). Bottom: distribution of hopping 
times; symbols correspond to simulation data at temperatures 
T = 1.0, 0.6, 0.4, 0.3, 0.25, 0.22, 0.2 (from left to right), and 
lines to fits using Eq. Q with fixed /3 = 1/2. 



a similar range of temperature and similar sizes. From 
Eq. (J), we have 



Tr(lnt) 



1 
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exp 



t 
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Symbols in the bottom panel of Fig. [3] indicate simula- 
tion data, while lines correspond to fits using Eq. |2"| 
with (3 — 1/2, i.e., the single fitting parameter is the re- 
laxation time i r ei(T). The fits are excellent over several 
orders of magnitude. Fig. 0] presents a similar analysis 
for the East model, for systems sizes L — 8£(T). The 
top panel shows C(t). The temperature dependence of 
the stretching is obvious, and follows well the theoreti- 
cal expectations. The bottom panel gives the data for 
hopping time distributions, and fits using Eq. I|2[). where 
now (3(T) is fixed from the corresponding values for C(t). 
Again, the fits are good over several orders of magnitude. 

The agreement between theory and simulations does 
not come as a surprise. This step was necessary, however, 
to establish the actual link between the dynamically het- 
erogeneous dynamics of the system and the distribution 
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FIG. 4: Same as Fig. 3 but for East model. Top: Temper- 
atures are T = 1.0, 0.6, 0.4, 0.3 (from left to right), and the 
stretching exponents are /3 = 0.55, 0.45, 0.37 and 0.27 (from 
left to right). Bottom: Temperatures are T = 1.0, 0.6, 0.4, 
0.3, 0.25 (from left to right), the fits are from Eq. @ with j3 
fixed by the correlators. 



of "happing" or "trapping" times as measured in Refs. 
[sL l9L fioj . It also shows that these names (hopping, trap- 
ping) somehow obscure the physical interpretation. 

IV. COMPARISON TO SIMULATIONS AND 
EXPERIMENTS 

The spin facilitated models studied above are a simple 
instance where the phenomenon of dynamical "trapping" 
occurs due to the presence of spatial dynamical hetero- 
geneities. Moreover, they are simple enough that the 
statistical properties of the "traps" can be worked out 
analytically. The basic initial observation of Stillinger 
and Weber was precisely the presence of these "hidden 
structures in liquid" 0, Il9| . More recently, quantitative 
studies of the statistical pro per ties of the IS dynamics in 
liquids have appeared H, U ^1 . The purely topographic 
interpretation of the data made in Refs. 0, M> HH was 
the main motivation for the current paper, since we think 
that the "basins" , "metabasins" or "traps" described in 
these works are a manifestation of the presence of dynam- 
ical heterogeneities, but not necessarily of static "traps" , 
as discussed in the previous sections. 



FIG. 5: Fits of LJ data from (top) and (bottom) 
using Eq. (|5J. The points are the published data, while the 
full lines are the fits, using two fitting parameters, f3 and t rc \ 
at each T. Top: (T,f3,t Tcl ) are (0.49, 0.225, 30641) (0.575, 
0.275, 272) (0.669, 0.316, 257) (0.764, 0.428, 135) from top to 
bottom. Bottom: (0.4, 0.234, 29500), (0.435, 0.254, 24200), 
(0.466, 0.246, 8820), (0.5, 0.210, 1268), (0.6, 0.358, 3138), 
(0.8, 0.380, 452), (1.0, 0.589, 203) from top to bottom. 



Having theoretical predictions from the spin facilitated 
models, it is tempting to reanalyze published data to 
check whether predictions are able to describe also more 
realistic models and experiments. Our ambition here is 
not to claim that one-dimensional spin models are enough 
to describe the physics of three-dimensional molecular su- 
percooled liquids, but rather to show that the mechanism 
described above is a simple and robust one, with an ap- 
plicability well beyond these simple models |2lj . 

In Fig. |SJ we present the data for hopping time dis- 
tributions from the simulations of Q (top) and |!j (bot- 
tom). The lines through the symbols are fits using the 
distribution for dynamical heterogeneities log times of 
the FA/East models, Eq. @. Apart from an irrelevant 
normalization constant, there are only two fitting param- 
eters, the stretching exponent (3 and the relaxation time 
t re i. The fits to the data of Denny et al. (top) are as good 
or better than those with a log normal distribution, as 
predicted by an activated scaling The agreement is 
also very good for Doliwa and Heuer's data (bottom), 
which we have chosen to present as 7r(lni) rather than 
as p(t) as they were originally published in |9(. The de- 
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FIG. 6: Distributions of molecular switching times in 
polyvinyl acetate at 299 K, from Ref. |23|. and fits using 
Eq. @ with p = 0.67, i re i = 0.070 s (Up), and /3 = 0.63, 
trei = 0.021 s (Down). 



viations in the fits in this case are a consequence of the 
data for the very short times investigated in |9j , where we 
expect the mechanism of dynamical facilitation to be less 
relevant, so that Eq. @ should describe the tails of the 
distributions more accurately than the small times |22j |. 

Finally, we have compared Eq. to experimental ob- 
servations. Although changing the system size is not 
conveniently realized in experiments, it is still possible 
to investigate the dynamics in "mesoscopic" regions, i.e., 
regions of a size comparable to the dynamic coherence 
length t(T). In Ref. [23|], time series for molecular po- 
larization were obtained by Vidal Russell and Israeloff in 
local probe experiments on polyvinyl acetate films. Typ- 
ical length scales probed through the use of non-contact 
atomic force microscopy techniques are ~ 20 nm, i.e., 
not very larg e compared to the estimated size of hetero- 
geneities Fortunately, since the polarization signal 
oscillates between two or four discrete levels, fast ther- 
mal fluctuations only produce oscillations about these 
discrete values, so that the corresponding distributions 
for switching times can directly be compared to Eq. J5J , 
therefore bypassing the IS construction. 

We reproduce in Fig. the distributions presented in 
the Fig. 5 of Ref. [23| , where we have adopted a log scale 
for the time axis, and represented 7r(lni) instead of pit), 
so that the comparison with the theoretical prediction 
@ can be made on the entire range of time scales and 
also to facilitate the comparison with the results pre- 
sented above. Again, we obtain an excellent fit of the 
experimental data to Eq. for the complete experi- 
mental time window. The stretching exponents we find, 
/3 = 0.67 and 0.63 for the up and down polarizations, 
res pec tively, are slightly larger than the ones estimated 
in [23. This is because the latter stemmed from a fit to 
a stretched exponential form of the tail of the distribu- 
tion. Instead, it is clear from Fig. [B] that the power law 
prefactor is necessary to account for the whole spectrum 



of time scales, and this shifts the stretching exponent to 
a slightly larger value. 

V. CONCLUSIONS 

We have shown that a consistent interpretation of the 
IS dynamics can be obtained solely invoking spatial and 
dynamical aspects related to the statistical properties of 
the dynamic heterogeneity of supercooled liquids. This 
physical picture was made more precise in the study of 
two simple dynamically facilitated spin models, which 
also allowed to obtain quantitative predictions which 
compare well to numerical and experimental data. Our 
results are additional evidence that dynamical hetero- 
geneities play a crucial role in the physics of liquids ap- 
proaching their glass transition [llj . 

In the topographic language, the so-called "landscape- 
influenced" regime evidenced by a decrease of (e/s) with 
decreasing temperature is the relevant one in numerical 
simulations [i. 12 1|. In our non-topographic perspective, 
this decrease of the IS energy corresponds to a decreas- 
ing concentration of facilitating defects |2l|, i.e., to the 
growth of the dynamical coherence length which is in 
turn directly responsible for the dynamical slowdown. It 
is therefore striking that the physics at relatively high 
temperatures, i.e., above the estimated mode-coupling 
temperature T c , can be described by the same physical 
mechanism of dynamic facilitation initially supposed to 
be relevant close to the experimental glass transition T g . 
Hence, our results imply that the dynamics above T c is 
heterogeneous, in agreement with simulations |25j. 

Moreover, as suggested in [3, it logically follows that 
the physics in this regime is qualitatively not very differ- 
ent from the one at temperatures closer to T g . This is 
at odds with the common belief that a change of mech- 
anism takes place close to T c This raises two inter- 
esting questions that we formulate as a conclusion to the 
paper. First, in the language of heterogeneities, what is 
the meaning, if any, of the crossovers reported close to 
T c ? Second, since these crossovers somehow justify the 
use of the mode-coupling theory above T c , why does the 
theory apparently work in a regime where the dynamics 
is heterogenous, and why does it eventually break down 
at lower temperatures? 
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